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The Anderson delocalization-localization transition is 
studied in multilayered systems with randomly placed inter- 
layer bonds of density p and strength t. In the absence of 
diagonal disorder (W = 0), following an appropriate pertur- 
bation expansion, we estimate the mean free paths in the 
main directions and verify by scaling of the conductance that 
the states remain extended for any finite p, despite the inter- 
layer disorder. In the presence of additional diagonal disorder 
(W > 0) we obtain an Anderson transition with critical dis- 
order W c and localization length exponent v independently 
of the direction. The critical conductance distribution P c (g) 
varies, however, for the parallel and the perpendicular direc- 
tions. The results are discussed in connection to disordered 
anisotropic materials. 



PACS numbers: 72.15.Rn, 71.30.+h, 74.25.Fy 



I. INTRODUCTION 

The understanding of the Anderson transition based 
on the scaling theory of localization jlj inspired many de- 
tailed numerical studies of disordered electronic systems 
The universality of the associated critical behavior 
was tested for various physical models, which include the 
crucial role of symmetry (^] and added magnetic field j| . 
Universal critical transport properties are also expected 
in the presence of hopping matrix elements which are 
not the same in the various lattice directions, as it can 
be seen from computations for weakly coupled chains and 
coupled planes ||[3j. It must be emphasized that many 
of the previous works on anisotropy include site diago- 
nal isotropic disorder and involve anisotropy only in the 
hopping magnitudes. This kind of anisotropy remains 
for zero disorder and is manifested in the band struc- 
ture. However, many realistic materials involve truly 
anisotropic disorder. For example, attempting to un- 
derstand the high-T c cuprates within a non-interacting 
electron picture in the presence of disorder requires the 
explanation of the contrasting resistivities in the paral- 
lel and perpendicular directions |^-|ic[|. Anisotropic site 
randomness in a form resembling a random superlattice 
with lateral inhomogeneities gave anisotropic localization 
for anisotropy below a critical value, even for arbitrarily 
small disorder ]Tl| . 

The in-plane resistivity for most of the layered 



high— T c materials exhibits metallic behavior, increas- 
ing linearly with temperature over a wide temperature 
range, while the perpendicular out-of-plane resistivity 
is very high at low temperatures and decreases rapidly 
as the temperature increases, reminiscent of semicon- 
ductors p2|-^5[. The contrasting behavior of the par- 
allel and the perpendicular resistivities was observed 
in P^S^-zLa^CuOy far below T c , down to the low- 
est experimental temperature |T(| . In the underdoped 
La2_ 3; Sr 3 ;Cu04 logarithmic divergencies of the corre- 
sponding resistivities accompanied by a nearly constant 
anisotropy ratio are, instead, observed suggesting an un- 
usual three-dimensional (3D) insulator. The electronic 
transport in these materials is expected to arise from 
scattering in the "insulating" layer between the conduct- 
ing CuC>2 layers [12 - 15 1 . On the other hand, in almost 



all high-T c cuprates doping impurities or oxygen vacan- 
cies occupy the insulating layers between the conducting 
"pure" Cu02 planes which implies interlayer disorder. 
Although the main aspects of transport in high-Tc ma- 
terials, such as the linear temperature dependence of the 
in-plane resistivity, are intimately connected with their 
strongly correlated nature, it is believed that anisotropic 
transport issues are, somehow, related to their layered 
structure. 

We propose a strongly anisotropic multilayered struc- 
ture (see Fig. 1) motivated by realistic anisotropic ma- 
terials. This system involves truly anisotropic interlayer 
disorder, anisotropic hoppings and the usual isotropic site 
diagonal disorder. It can be regarded as a very simple 
model for the cuprates where the CuC>2 planes are be- 
lieved to be identical without superlattice-like disorder. 
Our aim is to study both parallel (||) and perpendic- 
ular (_L) transport addressing the following main ques- 
tions: (1) does anisotropic localization occur (for exam- 
ple, localization in the layering direction and derealiza- 
tion within the layers) for interlayer disorder only? (2) 
with additional isotropic diagonal disorder is the critical 
behavior independent of the direction as scaling theory 
predicts? Firstly, we compute the conductance in the 
case of interlayer disorder alone to check whether its lo- 
calization behavior is the same in both directions. Sec- 
ondly, in the presence of additional isotropic disorder of 
strength W we obtain the critical disorder W c and the lo- 
calization length critical exponent v to see if they depend 
on the direction. We have also analyzed the statistical 
properties of the critical conductance distributions P c (g) 
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and although we can conjecture that it is a unique single- 
parameter function the distribution for the logarithm of 
the critical g in the "difficult" _L case resembles an insu- 
lator. 

In Section II we proceed with the definition of the 
tight-binding Hamiltonian for the multilayered lattice 
structure. In Section III, we consider the system with 
only interlayer disorder and estimate the corresponding 
mean free paths in the two directions. Our results in 
the absence of isotropic diagonal disorder allow to con- 
clude, in agreement with the scaling theory, that the 
states remain extended in both directions despite the 
strongly anisotropic interlayer disorder. However, the 
metallic conductance is very different for the _L case, be- 
ing insulator-like. In Section IV we review the numerical 
methods for the computation of the conductance in cubic 
and long wire systems. We find singular behavior along 
the layering direction due to the missing bonds. In order 
to avoid this problem we have developed appropriate nu- 
merical algorithms based on transfer matrix and Green 
function methods. Finally, in Section V we discuss the 
conclusions of the present study also in connection to 
realistic systems. 

II. RANDOM MULTILAYER LATTICE 

We propose a simple 3D anisotropic multilayered 
model which consists of parallel lattice planes randomly 
connected by interplane bonds as in Fig. I described by 
the Hamiltonian 

H = ^£ m ,i|m,Z >< m,/| + ^2 (|m,Zxm',Z| 

ml (m,m'),l 

+t' ml \m,l >< m,l + 1| +H. a), (1) 

where m, m' denote the two-dimensional site indices in 
each layer and I is the layer index. The first term in 
Eq. (I) describes diagonal (isotropic) disorder with the 
site matrix elements e chosen randomly from a box dis- 
tribution within [-W/2, W/2], the second term describes 
nearest-neighbor hopping of unit strength within the lay- 
ers, which sets the energy unit, and the third term corre- 
sponds to interplane hoppings t' ml = or t, placed with 
probability p at random layer positions m. The inter- 
plane term obeys the binary distribution 

P(C,i) = P*(C,i,*) + (1 - PWC,i,0). (2) 

The proposed structure has both anisotropic hoppings 
due to t and anisotropic disorder due to p. The trans- 
port characteristics are obtained by calculating the con- 
ductance along the || and the _L directions. 

The missing perpendicular bonds in the layering di- 
rection disturb particle migration even the presence of 



interlayer disorder alone. Unlike a naive expectation we 
find no critical point when W — 0, for any p ^ 0. We 
show that the system is metallic independently of the 
direction, although the behavior of the conductance is 
very different in the two directions. In the presence of 
additional diagonal disorder, denoted by W, a critical 
disorder W c is obtained for various choices of the density 
p and strength t. The critical point W c within finite size 
errors is found to be the same in both directions. 

III. MEAN FREE PATHS FOR INTERLAYER 
DISORDER 

We consider anisotropic disorder in the perpendicu- 
lar layering direction represented by p and t, due to 
the randomly placed bonds among consecutive layers, 
in the absence of diagonal disorder W. In this disor- 
dered anisotropic lattice one might expect transport to 
be hindered in the perpendicular direction. It is worth 
examining whether is present or not. In order to pro- 
ceed we adopt a convenient layer-diagonal representation 
since for W — the 2D layers are perfect planes and can 
be easily diagonalised. The eigenstates at the I— th layer 
|k||,Z) are labelled by the parallel momentum k|| and the 
Hamiltonian H can be expressed in the convenient Bloch- 
Wannier basis 

|k||,0 = -^=]>> 4k ''- m |m,0, (3) 

with parallel momentum k|| , the layer index I and m 
summed over all N\\ — L 2 sites in every layer for a systen 
with L 3 sites. For p = the 2D layers are perfect and 
k|| is a good quantum number. For p = 1 the system 
reduces to a perfect 3D lattice and both k|| , k± become 
good quantum numbers. 

We consider the case of p ^ 0, 1 where the translational 
symmetry in the plane directions is also broken and k|| 
is no longer a good quantum number. In this mixed 
representation the Hamiltonian H can be expressed as 

F = ]Te||(k||)|k||,0<k||,Z| 
k,|,i 

+ E E^ml k lh z )( k ll' Z + 1 l +H - c -]' ( 4 ) 
( k ii. k y 

with the parallel kinetic energy 

e||(k||) = 2 cos(fc x ) + 2cos(fcj / ) (5) 

and the hopping matrix element between neighboring 
planes 

<m=iv-E el(kll " k|l) - mi U- (6) 
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In order to investigate the question of localization in the 
layering direction we define the retarded Green function 

GO^ljkf,,*';*) = -te(t)<[c k „ ) «(t),4, |iI ,(0)]+) J (7) 

where Ck||,z(i) is the time-dependent destruction opera- 
tor of electron in the state |k||,/). Its diagonal element 
G(k|| , I; k|| , I; t) gives the probability for finding an elec- 
tron on the layer I with momentum k|| at time t, if ini- 
tially it was on the same layer having the same momen- 
tum. The Fourier transformation of the diagonal Green 
function with respect to time is 



G(ky,I;k,|,I;E) 



1 



E - E(kii , I, E) 



with the self-energy Sii written as |lj 



Z(k h l,E)=e ll (k ll )+J2J2 T ; 

71=1 j 



(n) 



(8) 



(9) 



The ro-th order term is a sum over all paths j of length 
n starting and ending in the same state |k|| , 1} 



if>=t Ml (k lhk||1 )n 



.k||i,k|| 



t 



(10) 



tti E -H^\\i) 
with k||„ +1 = k|| and l n +i = I- From Eq. (|^) one has 



meBi. 



with the sum for m over the set of sites Bi u i i+l which have 
interplane bond connections. It is seen from Eq. JTl] ) 
that the diagonal in momentum matrix elements (kiu = 
k|| i+1 ) are exactly pt. For a given order n if E lies within 
the pure 2D band (k|| ) the most divergent term of Eq. 
(^) comes from the path which has intermediate states 
k||j,Zi) with en (kiu) ~ E for all i. The corresponding 
term approaches {pt/[E — e^k^)]}™ which sequentially 
connects nearest neighbor plane states with the same kiu. 
It can be also seen that momentum scattering is always 
accompanied by interlayer hopping since a change of k|| 
leads to a change of the layer index. 

The above analysis splits the Hamiltonian into two 
parts H = H + Hi . The "undisturbed" part H repre- 
sents a perfect anisotropic 31? lattice with intralayer (in- 
terlayer) hopping 1 (pt) . The rest is a "random" Hamilto- 
nian H\ with only off-diagonal matrix elements expressed 
in the iJo-diagonal basis |k) = |k||, k z ) of the form 



H!(k,k') = — '- 

N 

for k|| ^k'y, 

#l(k,k') = 0, for k,, = k',, 



(ky — kj, )-m+i(k z -h' z )l 



(12) 



where N is the total number of lattice sites and k z the 
perpendicular momentum. 

It can be seen from Eq. ( p^ ) that the matrix elements 
of Hi are complex numbers of average amplitude value 
modulo 1 plus a random phase. If the size of the system 
increases to infinity the phase exhausts all possible values 
in [0, 2ir] and the average should vanish. In this situation 
the scattering by the random configurations of the inter- 
plane bonds can be well described by perturbation theory 
with k-space self-energy 



£(k,£)^ e (k)+]T 



|gi(k,k')| 2 
E - e(k') - 



(13) 



The configuration average (Im^ k , E^}nj)-ir ) can ^ e 
computed as a function of p and the results fitted to a 
semicircular form as p(E)t 2 p(\—p), with E within the Hq 
band and p(E) the corresponding density of states. This 
allows to estimate the lifetime of states r ~ — r 
and the corresponding mean free paths 



An = 



1 



p(E)t 2 p(l-p) v / 2+ P H^ 



p(E)(l-pW2+pH2 



(14) 



(15) 



for Fermi velocites 



IK-L)- 



We observe that for small-p the obtained mean free 
path in the parallel (perpendicular) direction is pro- 
portional to 1/p (p). This implies that the scattering 
strength of the interlayer disorder increases with p for 
parallel transport but decreases with p for perpendicular 
transport. This intrinsically anisotropic situation seems 
in contradiction with the scaling theory of localization 
because p — * might be thought to be in favor of local- 
ization in the perpendicular direction, while transport is 
the least affected giving extended states in the parallel di- 
rection. However, the application of the one-parameter 
scaling theory to this situation should give a common 
critical point in all directions so that transport in the 
perpendicular direction should be extended as well. This 
is, indeed, numerically confirmed in Chapter IV where 
we show that the scaling behavior of the conductance in 
finite cubic systems with interlayer disorder is in agree- 
ment with the scaling theory. It turns out that the in- 
terlayer disorder is not sufficient to localize the electrons, 
even in the "difficult" layering direction, without any ad- 
ditional diagonal disorder W. The statistical behavior 
of the conductance, however, is very different in the two 
directions. 
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IV. NUMERICAL CALCULATION OF THE 
CONDUCTANCE 

A. The cube 

The parallel (perpendicular) conductance ^r9\\(±)(L) 
can be obtained foraLxIxL cubic system at the Fermi 
energy E directly from the multichannel Landauer- 
Buttiker formula Jlf] 

ff || (±3 (i)=Tr(t+ ±) t|| (±) ), (16) 

where tiinj is the transmission matrix for electronic 
propagation along the || (_L) direction computed by trans- 
fer matrix techniques. Two perfect semi-infinite bars are 
attached to two opposite sides of this cube and hard-wall 
boundary conditions are used for the other sides. The 
number of independent channels for the incoming and 
outgoing leads is L 2 and the transmission matrix can be 
calculated from the amplitudes of transmitted waves in 
the outgoing leads by assigning, at a time, a unit incident 
wave amplitude for one channel in the incoming leads and 
zero for the rest. We can easily establish the recursion 
relations for the corresponding wave function coefficients 
along the parallel direction to calculate g\\. In the _L 
direction computing the matrix tj_ is not possible since 
the recursion relations are singular due to the presence 
of zero hoppings for the missing interplane bonds. To 
overcome this difficulty we set up the recursion relations 
along the direction parallel to the planes, perpendicular 
to the leads, but with the reflection and transmission 
coefficients in the channels as the unknown variables. 
We can solve these recursion relations with hard wall 
boundary conditions perpendicular to the leads. By this 
method we obtain the perpendicular transmission matrix 
t_L avoiding the singularities due to missing bonds. This 
is a convenient tool to consider propagation in the layer- 
ing _L direction by a transfer matrix product only along 
the easy || direction. In order to suppress fluctuations 
we have taken averages over up to a 5000 random cubic 
configurations in each case JToh ■ 

B. The wire 

We can also compute the parallel (perpendicular) 
dimensionless conductance g\\ (M) for a quasi-one- 
dimensional M x M X L geometry, via Green function 
methods [p0[ . In the parallel direction m = x,y the 
Hamiltonian of the M x M slice is incorporated into the 
transfer matrix 

rp _(V-l +1 Q\(E-H x -l\(l \ 
x ~ \ l)\ 1 o ){o V x - hx )> 

(17) 



where the matrix V x , x +i has unit elements. For large 
length L in the x-direction the product 

L 

T = [] T x , (18) 

x =i 

has eigenvalues exp7^L and Lyapunov exponents 7,, i = 
1,2,..., M 2 . The smallest positive Lyapunov exponent 
71 determines the scaling parameter A||, via 

V = ^ = (m)" (19) 

with £m = the largest localization length, which is of 
interest for finite-size scaling studies when the width M 2 
of the slice is varied. We obtain the critical value A c || 
at the point where A11 becomes independent of M (see 
Table 1). In order to ensure accuracy of about 1% for 
A || the length of the studied system is more than about 
200000. We find that this length near the critical point 
varies as oc A -1 . In the case of very strong anisotropy we 
were unable to obtain this accuracy for all M and W. 

This method cannot be used in the ^-direction. The 
reason is, again, the zero elements of the hopping matrix 
Vy+i with probability 1 — p so that the inverse matrix 
^ll+i w hich enters ( |l7| ) becomes singular. In order to 
avoid this problem we can alternatively use the Green 
function G(E) = (E — H)^ 1 by applying the iterative 
scheme of J2(| based on the two equations 

r (l + l) _ r {l) v r (l + l) / 9r) N 

W,z+i - u 1,1^,1+1^1+1,1+1 V Z[} ) 

and 

G^ + i = \E- m+i - V^uGfjVu+i]- 1 - (21) 

The Hamiltonian Hi represents the Zth M x M layer, 
Vi,i+i the hopping between layers 1,1 + 1 and Gf \ 2 is 
the Green function of the system with length I between 
layers l\ and li- The diagonal matrix Vi t i+\ of the order 
of M 2 has zeros with probability 1 — p so that V{~ l+X 
becomes singular. The main advantage of the formulae 
( pp| , ^l|) is that they do not contain the inverse of Vi,i+\. 
The disadvantage is the necessity to invert a matrix of 
order M 2 [Eq. (pi])] at each iteration step. In this case 
we restrict the number of iterations to L < 40000 and the 
corresponding peprendicular scaling parameter is defined 
as 

Al^f logTrGf L = (^) _1 . (22) 

The critical disorder W c and the critical exponent v are 
computed in both directions from the numerical data of 
A(M, W). The usual linearization near the critical point 

log A(M, W) = a M + /3m log W, (23) 
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and the independence of A on M at the critical point 
gives W c from the slope of the linear dependence 



a.M = - log W c x (5m + const. 



(24) 



The obtained critical A c y and A c _l are different from the 
value A c of the corresponding isotropic system (see how- 
ever Eq. (|2S| ) below). The critical exponent v is deter- 
mined from the coefficients /3m via 



log Af 

"log/3 



M 



(25) 



We have also repeated the computations for cubic sys- 
tems and although the finite size effects became more 
pronounced the obtained critical values differ very little 
from those obtained for long wires. 

In agreement with the W — case it is reasonable to 
suppose that the form of P c (g) also depends on the di- 
rection. For the || direction the conductance is calculated 
from the formula 



g(M) = £ 



^cosh^/2) 



(26) 



where the Zj's are the logarithms of the ith eigenvalues 
of the matrix t T t 0] and in the limit L » M converge 
to 2Mji. In the _L direction the singular behavior of the 
matrix Vy+i does not permit to use the formula (^6|). To 
calculate log g± in this case we use the fact that A c ^ < < 
1 for most of the critical points discussed below. Then, 
the critical conductance from z\ = 2/A c is very small 
and can be estimated from the contribution of the first 
channel via 



logg± = Tr G hL (E + iri)G L ^{E - it)), 



(27) 



with the imaginary part of the energy r\ given by the ratio 
of the bandwidth over the mean level spacing E2] . 



V. RESULTS 

A. W = 

Fig. 2 shows the scaling behavior of g\\r±)(L) for the 
L 3 -site cubic system for various bond densities p and 
anisotropic interplane coupling t = 0.3, in the absence of 
diagonal disorder W . The parallel conductance is shown 
to increase ballistically (~ L 2 ) for small L and linearly 
for higher L. In the large size (L — > oo) limit the cor- 
responding scaling function 0(L) — d\ng/dlnL becomes 
positive for <?||(_l), which implies extended states in both 
directions for any p, in agreement with the scaling the- 
ory which predicts a common critical point in any direc- 
tion. However, the obtained transport behavior is essen- 
tially different in the two directions. In Fig. 3 we show 



the energy dependence of the conductance for p = 0.5 
and t — 0.3 where the ratio of the two conductances is 
close to the estimate g\\/g± ~ (t\\/t±) 2 for t\\ = 1 and 
t± = pt = 0.15 Q. A key finding from Fig. 2 is a rather 
smooth 0ii (E) while g±(E) displays violent oscillations as 
a function of E. The dips in g±(E) can be regarded as 
due to "minigaps" in the perpendicular direction which 
might have effect similar to a semiconductor, leading to 
insulating kind of behavior for the out-of-plane conduc- 
tivity when the Fermi energy is varied. 



B. W > 

The W^-dependence of A(M, W) for different M's and 
various parameters p, t is presented in Figs 4., 5., 6., 7. 
The corresponding critical points are calculated by the 
described fitting procedure. The results are listed in Ta- 
ble 1 and the data are very reliable for the || direction. 
In the _L direction they arc much harder to analyze due 
to finite-size effects. For example, in case D (Fig. 7) the 
obtained M-dependence of A±(M) is not monotonic and 
for W = 5 we find A(M) which decreases with M for 
small system sizes, imitating insulating behavior. How- 
ever, for larger M > 12 the values of A(M) begin to 
increase and the correct scaling is restored. This in- 
dependence is caused by a second irrelevant scaling pa- 
rameter from the relation A = aM 1/l/ + bM 13 with j3 < 0. 
The correct estimation of the critical parameters in the _L 
direction requires either numerical data for larger M or 
possibly more sophisticated fits |23| . Although for strong 
anisotropy we could not obtain accurately the critical pa- 
rameters in the _L direction we check that our results 
converge to those of the || direction when M grows. 





P 


t 


W c 


A c 


V 


M 


(5) 


var g 


1 




A 


0.6 


1.0 


14.47 ± 1.10 


0.676 


1.619 ± 0.14 


10-16 


0.59 


0.16 


B 


0.6 


0.3 


10.48 ± 0.58 


0.933 


1.666 ± 0.10 


6-14 


1.5 


0.4 


C 


0.1 


0.3 


7.93 ± 0.24 


1.218 


1.514 ± 0.06 


12-18 


2.6 


0.7 


D 


0.6 


0.1 


8.05 ± 0.45 


1.174 


1.790 ± 0.14 


6-16 


2.6 


0.7 


1 


A 


0.6 


1.0 


14.30 ± 1.30 


0.437 


1.546 ± 0.16 


6-14 






B 


0.6 


0.3 


10.20 ± 0.67 


0.189 


1.308 ± 0.21 


8-18 






C 


0.1 


0.3 


7.18 ± 0.74 


0.089 




12-20 






D 


0.6 


0.1 


6.80 ±1.94 


0.100 




14-18 







TABLE I. Review of various critical parameters with the 
mean conductance (g) and its variance var g = (g 2 ) — (g) 2 . 
The obtained values confirm that the critical disorder is 
very close for || and _L transport, in agreement with the 
one-parameter scaling theory. We could not calculate the crit- 
ical exponent v in the _L direction for the cases C, D. 



5 



In this direction the critical region is very narrow and 
we could neither calculate the critical exponent, since 
larger system sizes are required. Nevertheless, the scaling 
analysis for A and B gives satisfactory results in both 
directions which confirm the equal W c and v. Moreover, 
we find the quantity 

A c = [A*,, x A C J 1/3 (28) 

which gives the critical value of the corresponding 
isotropic model H. 

Fig. 8(a), (b) presents the probability distribution of 
the critical conductance P c (<7|i) in the parallel direction. 
The distribution is shown to be size-independent but de- 
pends on the various critical points. For C and D the 
obtained P c (<7|i) and A c || are the same. This indicates 
that P c (<7||) is determined only by A c || |^|. The numer- 
ical data for the mean and variance of the conductance 
are also unique functions of A c m (see Table 1) supporting 
this conjecture. In the limit t — > the critical distribu- 
tion in the parallel direction is expected to converge to a 
Gaussian. However, the spectrum of the obtained higher 
Lyapunov exponents shows square-root behavior similar 
to isotropic 3D disordered systems p2|| . To display the 
dramatic differences in parallel and perpendicular trans- 
port, we also present the critical distribution of logg±, 
calculated for the perpendicular direction for the criti- 
cal case C. This distribution has all the features of the 
localized regime since it is log-normal with 

var log g x ~ - (log g±). (29) 

The important difference with the insulating regime is 
the fact that P c (loggjJ remains system-size invariant. 
The relation of Eq. (|2^) is also valid for the critical points 
B and D where in the _L direction (g±) << 1. 

VI. DISCUSSION 

The random topological multilayered structure studied 
may be regarded as a first step towards an explanation, 
via non-interacting electrons, for properties of strongly 
anisotropic materials. For example, in the case of W = 
the 2D layers are perfect and the disorder represented 
by p can be due to impurities or oxygen vacancies in the 
insulating layer among the 2D planes of the cuprates. 
The electrons propagating in the perpendicular direction 
of this system shall encounter anisotropy in the disor- 
der due to the random interplane links in addition to 
the value of t which can be different to that of the par- 
allel direction. However, the distribution of the critical 
conductance P c {g) depends both on the choice of the pa- 
rameters and the direction where the electron moves. 

The considered anisotropic structure exhibits rather 
strange transport properties on a given scale, expressed 



in the dramatic differences of the critical conductance in 
the parallel and perpendicular directions. In the perpen- 
dicular direction the conductance distribution slightly 
bellow the critical point is log-normal resembling the sta- 
tistical properties typical of an insulator. A similar sta- 
tistical "anomaly" has been described in [^0|. However, 
a strong difference to a "true" insulating regime exists 
since the conductance still grows with the size. In the 
large size limit the corresponding distribution reaches a 
Gaussian. An analogous discussion holds for the paral- 
lel direction slightly above the critical point. The main 
criterion for the specification of the critical regime is the 
size dependence of the conductance (g) (or (\ogg)). 

The proposed model for W = may have some re- 
lation to the strongly anisotropic transport proprties 
of high-T c cuprates. In these materials, as tempera- 
ture increases, the inelastic scattering due to phonons, 
spin waves or other excitations within the CuO planes, 
can cause a decrease of the inelastic scattering length 
li n . If the temperature is so high that 2j„ becomes 
smaller than mean free path the transverse conductiv- 
ity is metallic. Experiments for Bi2Sr2CaCu20s and 
underdoped I^-xSr^CuO^ YBa2Cu306+;c give out-of- 
plane resistivity which has a semiconductor-like temper- 
ature dependence at low temperatures (high at small- 
T with a rapid decrease by increasing T) and a linear- 
in— T behavior at high temperatures. The characteristic 
crossover temperature between the two regimes T* de- 
creases by increasing the doping in La2- a; Sr a ;Cu04 and 
YBa2Cu306+a; ]7|-[l0|]. We notice that if we relate the 
bond density p with the doping density of the high-T c 
materials the obtained p-dependence of the perpendicu- 
lar mean free path can be used to explain qualitatively 
the reported behavior. It must be pointed out that the 
relation between the cuprate doping density and the bond 
density p is natural, since an increase in the number of 
the doping impurities or the oxygen atoms in the layer 
between two Cu02 planes increases the number of hop- 
ping paths between the two planes. As T* decreases fur- 
ther (below T c ) the out-of-plane normal-state resistivity 
also becomes metallic, which has been observed in high- 
quality single crystals of YBA2CU3O7 and other high-T c 
cuprates corresponding to the absence of disorder (p«l 
with W — in the proposed model) with almost infinite 
perpendicular mean free path. This effect occurs only in 
the perpendicular direction since the parallel mean free 
path is always much longer and inelasting scattering be- 
comes dominant. 

In summary, we have introduced a simple layered lat- 
tice model with anisotropic disorder described by the 
interplane bond density p, in addition to the usual 
anisotropic band structure expressed via the interplane 
hopping t. In the absence of diagonal disorder we show 
extended states in both directions but the obtained mean 
free path and the conductance in the _L direction is much 
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smaller than in the || direction. Moreover, g± fluctu- 
ates strongly as a function of energy, which leads to an 
insulator-like temperature dependence of the conductiv- 
ity in the _L direction. In the presence of additional di- 
agonal disorder of strength W we have shown that the 
critical disorder and the critical exponent v do not de- 
pend on the transport direction. The obtained data for 
the localization exponent v agree with recent accurate 
estimates for the isotropic model p5| and confirm the 
universality at the metal-insulator transition. The ob- 
tained critical conductance distribution P c {g) although 
independent on the system size depends strongly on the 
parameters and the direction of transport. 
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FIG. 1. A picture of the multilayered structure which consists 
of 2D square lattices I (layers) connected by perpendicular bonds 
of strength t placed at random positions with probability p. 



FIG. 3. (a) The energy-dependent g\\ for a cubic layered system 
with L=10,15 and W = 0, t=0.3, p=0.5. (b) The g± is much 
smaller and exhibits violent oscillations as a function of the energy 
E. 
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FIG. 2. The gy as a function of the linear system size L for a 
cubic Lx Lx L system of parallel planes with W = and randomly 
placed interplane bonds of density p with strength 4=0.3. In the 
inset g± for the same system exhibits similar behavior but much 
smaller values. 
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FIG. 4. A: The behavior of the scaled localization length 
S.M /M for the parallel and perpendicular direction in the M X M xL 
system with p = 0.6 and t = 0.1 as a function of W. The critical 
point is located at W c ~ 14.47 in the || direction and W c ^ 14.30 
in the X direction. 
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FIG. 6. C: The behavior of the scaled localization length §m/A^ 
for the parallel and perpendicular direction in the MxMxL system 
for p = 0.1 and t = 0.3 as a function of W. The critical point is 
displayed at W c ~ 7.93 in the || and W c ~ 7.18 in the _L direction. 
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FIG. 5. B: The behavior of the scaled localization length §m/.M 
for the parallel and perpendicular direction in the M X M X L system 
with p = 0.6 and t = 0.3 as a function of W. The critical point is 
located at W c ~ 10.48 in the II and W c ~ 10.20 in the _L direction. 
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FIG. 7. D: The behavior of the scaled localization length 
£,m/M for the || and _L direction in the M X M X L system with 
p = 0.6 and t = 0.1 as a function of W. The critical point is located 
in Wc ~ 8.05 in the || and Wc ~ 6.80 in the _L direction. It is seen 
that the data for smaller M = 6,8 fail to cross at the same point 
indicating "insulating" behavior. 
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Log(g 1 )/<Log(g_ L )> 

FIG. 8. (a). The critical distribution of the parallel conduc- 
tance in the case A with < grii >R3 0.59, var(gy) = 0.16. (b). The 
same as in (a) for the critical points C (full symbols) and D(opcn 
symbols) with < gn >fti 2.6, var (gn) = 0.7. The C,D have the 
same A c and the same critical distribution (see Table 1). (c). The 
critical distribution of P(loggx) m the perpendicular direction for 
case C is also shown for comparison. 
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